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ABSTRACT 

The flame born in the deep interior of a white dwarf that becomes a Type la supernova 
is subject to several instabilities, the combination of which determines the observational 
characteristics of the explosion. We briefly review these instabilities and discuss the length 
scales for which each dominates. Their cumulative effect is to accelerate the speed of the 
flame beyond its laminar value, but that acceleration has uncertain time and angle dependence 
which has allowed numerous solutions to be proposed (e.g., deflagration, delayed detonation, 
pulsational deflagration, and pulsational detonation). We discuss the conditions necessary for 
each of these events and the attendant uncertainties. A grid of critical masses for detonation 
in the range fO^ - 2 x fO^ g cm~'^ is calculated and its sensitivity to composition explored. 
The conditions for prompt detonation are discussed. Such explosions are physically improbable 
and appear unlikely on observational grounds. Simple deflagrations require some means of 
boosting the flame speed beyond what currently exists in the literature. "Active turbulent 
combustion" and multi-point ignition are presented as two plausible ways of doing this. A 
deflagration that moves at the "Sharp- Wheeler" speed, O.lgcst, is calculated in one dimension 
and shows that a healthy explosion is possible in a simple deflagration if fuel can be efficiently 
burned behind a front that moves with the speed of the fastest floating bubbles generated by 
the non-linear Rayleigh- Taylor instability. The relevance of the transition to the "distributed 
regime" of turbulent nuclear burning is discussed for delayed and pulsational detonations. This 
happens when the flame speed has slowed to the point that turbulence can actually penetrate 
the flame thickness and may be advantageous for producing the high fuel temperatures and 
gentle temperature gradients necessary for detonation. No model emerges without difficulties, 
but detonation in the distributed regime is plausible, will produce intermediate mass elements, 
and warrants further study. The other two leading models, simple deflagration and pulsational 
detonation, are mutually exclusive. 



Subject headings: hydrodynamics - stars: supernovae: general - turbulence 



1. Introduction 

Despite 25 years of intensive investigation (e.g., Arnett f969), the basic physics whereby a carbon- 
oxygen core of nearly the Chandrasekhar Mass (f .39 Mq) explodes as a Type la supernova (SN la) is still 
debated. One may reasonably conclude that it is a hard problem. In fact, only recently has the astrophysics 
community begun to profit from the extensive experience of the chemical combustion community in order 
to appreciate fully just how complicated burning coupled to hydrodynamics really can be. 
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The astrophysical problem is especially hard because the nuclear flame propagates in an extensive 
medium in which gravity plays a role and several instabilities have time to develop over a large range 
of length scales. Any realistic solution must take cognizance of these instabilities and, if only by 
parameterization, incorporate them into the stellar model. 

This need, and variability in the outcome depending on uncertain parameters, has given rise to several 
classes of supernova models, all largely empirical. These include the "delayed detonation" (Khokhlov 
1991abc; Woosley & Weaver 1994), "pulsational detonation" (Arnett & Livne 1994ab); "pulsational 
deflagration" (Nomoto, Sugimoto, & Neo 1976; Ivanova, Imshennik, & Chechetkin 1974); "convective 
deflagration" (Nomoto, Thielemann, & Yokoi 1984); and the fractal model of Woosley (1990), each of 
which, by various contrivances, generates a flame which is born slow and accelerates very rapidly as the star 
begins to come apart. This behavior has been found essential to obtaining nucleosynthesis that agrees with 
spectroscopic observations. 

There are two simple solutions to the explosion problem, neither of which is thought to be correct, 
but which bound the real solution - a) a laminar conductive flame, and b) prompt detonation. The latter 
is improbable ( 3.2.2| ) and would give unacceptable nucleosynthesis - no intermediate mass elements; the 



former is unphysical and would not give an energetic explosion. Between these two extremes lies the 
regime of unstable flame propagation to which this paper is mostly devoted. We begin by reviewing 
the relevant instabilities and current views regarding their importance and mutual couplings. Though 
much has been written, especially regarding the Rayleigh- Taylor (RT) instability, not all share the same 
views even on this fundamental subject. It is therefore necessary to state (our view of) the basics before 
proceeding to the models, and so we briefly discuss (||) the RT-instability, the Landau-Darrieus (LD) 
instability, the Kelvin-Helmholtz (KH) instability, and turbulence. Our main goals here are: 1) ascertaining 
the uncertainties in the models in an attempt to resolve the leading candidate(s); and 2) exploring the 
implications of new physics - active turbulent combustion, multipoint ignition, and, especially, detonation 
in the distributed regime - for models for SN la's. 

It is frequently stated (Khokhlov 1995, Arnett & Livne 1994a) that simple deflagrations, those in which 
the flame remains at all times subsonic and in which an intervening pulse of the white dwarf does not occur, 
cannot give an energetic explosion. We discuss why this may not be true, either because of physics that 



has been left out of the calculations ( 3.1.2 , 3.1.3 ) or because the flame has not been modeled sufficiently 



accurately in three dimensions (3.1.1) 



Another important point that has not been adequately discussed in the context of astrophysical flames, 
is the idea of "distributed burning" . So long as the thickness of the flame is negligible, it can be treated 
as a temperature discontinuity. Heat never moves far from the burning surface. This is the case during 
most of the flame's life. But for densities below ^ 3 x 10^ g cm~^ turbulence will disrupt the flame sheet. 
Once this happens, one can no longer speak of the flame as either laminar or conductive. It has entered 
a domain where turbulent transport directly affects the nuclear reaction region and smears it out over 
macroscopic length scales. The combustion community calls this type of burning the "distributed flame 
regime" (e.g., Pope 1987) and it is a frontier topic for them as well as us. While the increased mixing of 
hot ash and cold fuel can, in some instances, be beneficial for provoking a detonation, the same process also 



leads to compositional mixing which increases the critical mass required for detonation to occur (3.2.1). In 



section (3.2.5) we discuss the physical conditions required for this sort of detonation and show that it is a 



reasonable, if uncertain, occurrence. In this sense, it is superior to many other models in the literature. 

Much has been written about other forms of the "delayed detonation model" (e.g., Khokhlov 1991abc; 
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Woosley & Weaver 1994), yet the physics of the transition to detonation remains obscure. In sections (3.2.4) 
and ( ^.2.5| ), we discuss the two different kinds of delayed detonation that have been previously published 
and why each is unlikely to occur. 



In the conclusions l4h we summarize our results. 



2. FLAME DYNAMICS AT HIGH DENSITIES 

2.1. The conductive laminar flame 

The simplest solution to the propagation of burning in a premixed fuel is an elementary "flame" . Heat 
is transported ahead of the burning region, in this case, by electron conduction. The temperature rises 
to the point where reactions can consume the fuel (carbon) on a diffusive time scale and this condition 
sets both the thickness of the flame and its steady velocity (Landau & Lifshitz 1991). In the absence 
of instabilities, these quantities can be determined analytically with considerable precision (Timmes & 
Woosley 1992). 

The critical mass required to keep a flame alive is small, a few times 4/37rpZj^, with Zth, the flame 
thickness. For isolated regions below this mass, heat can diffuse out and the flame will die, but so long as 
a critical mass remains intact, one cannot extract heat from the ashes of the combustion (in order to raise 
the temperature of the fuel) over a greater distance than /th- It is impossible for a simple laminar flame to 
turn into a detonation. So long as there is a flame with a well defined surface (exceptions to this will be 
discussed in sections |2.6| and 3.2.5| ), detonation can only be achieved by increasing the area of the flame. 



It is well known that the laminar speed is too slow to make a supernova with the observed properties. 
In fact, we shall conclude (see also Niemeyer & Hillebrandt 1995b and Khokhlov 1995) that the laminar 
speed is not very relevant for the effective rate at which burning spreads. Turbulence and Rayleigh- Taylor 
instabilities are more important and carry the flame at a speed independent of the microphysics. However, 
it is important to keep track of the flame thickness, as this may ultimately affect even the macroscopic 
nature of the burning, and to distinguish the speed with which the outer boundary of burning spreads 
from the total rate of mass consumption. These rates are quite different if a large amount of fuel becomes 
entrained. 

After a brief review of the various instabilities that arise from a linear analysis of flame propagation 



under the influence of gravity and shear in section (2.2), we discuss the quasi-stationary structure of the 



flame surface after it reaches the fully nonlinear regime in sections (p.3[), (2.4) and (k 



2.2. Linear instabilities 

The problem of linear hydrodynamical stability of subsonic flames in the thin flame representation 
was first analyzed by Darrieus (1938) and, independently, by Landau (1944). Landau's result for the linear 
growth rate also includes the influence of gravity on the density jump produced by the flame front, which 
is equivalent to the Rayleigh- Taylor instability (Chandrasekhar 1961) if the flame speed is set to zero. In 
the context of nuclear flames in degenerate carbon, the LD instability has been explored in detail both 
analytically and numerically by Blinnikov & Sasorov (1996) and Bychkov & Liberman (1995), and its onset 
has been demonstrated by means of two dimensional hydrodynamical simulations (Niemeyer & Hillebrandt 
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1995a). All these studies agree that on scales larger than the Markstein, or critical, length ^crit ~ 100 Zth 
(Markstein 1951), flames moving upward are unstable to both LD and RT instabilities on all wavelengths. 

In addition to radial perturbations, we have to consider perturbations that face in a direction 
perpendicular to gravity, where buoyancy of burned material floating in the cold fuel induces tangential 
velocity differences along the flame surface. Here, another instability becomes important, the KH or shear 
instability (e.g.. Landau & Lifshitz 1991). In our context, it is quite important to know the circumstances 
under which the flame surface in a Type la supernova is KH unstable since this condition marks the 
transition to the fully turbulent burning regime. Strictly speaking, tangential discontinuities can only occur 
if the mass flux through the surface vanishes, for otherwise continuity of the momentum flux imposes 
continuity of the tangential velocity components. The mass flux through the burning zone is by definition 
non-zero. Thus, we can qualitatively argue that the propagation of a flame tends to stabilize the front 
against KH-instability. If, however, the flow field around a burning bubble is dominated by buoyant 
acceleration, i.e. the mass flux becomes small compared with the velocity components tangential to the 
front, the flame becomes KH unstable, as shown by numerical simulations (Niemeyer & Hillebrandt 1996). 
This occurs during the nonlinear stage of the RT-instability. 



2.3. The cellular regime 

In the nonlinear evolution of the LD instability, the formation of cusps at the points where the flame 
front self-intersects gives rise to an additional quadratic damping term for the perturbation amplitude that 
is not included in Landau's linear stability analysis (Zeldovich et al. 1985). It stabilizes the flame surface 
after cells with a stationary, scale independent amplitude have formed. The speed of the emerging cellular 
surface is directly given by the increased surface area, yielding 

Mccll = Mlam(l + 4p)) , (1) 

where the velocity increment e is a function of the gas expansion parameter fj, = py^/p^ (Zeldovich 1966). If 
H is small, as in the case of burning in degenerate matter, e can be approximated by (Zeldovich 1966) 

= 0.41(1 -M)^ (2) 

Motivated by the large dynamical range of thermonuclear flames in white dwarfs Blinnikov and Sasorov 
(1996) proposed a fractal model for the cellular structure of LD unstable flames. By means of a statistical 
analysis of the Sivashinsky equation for thin flame propagation, the authors estimated the fractal dimension 
of one dimensional flames as 

Did = 1 + Doj^ (3) 

where 7 = (1 — /u). This result was confirmed by numerical simulations of the closely related Frankel equation, 
which yielded Dq rs 0.3. Furthermore, the authors derived the fractal dimension of two-dimensional fiame 
surfaces as ~ 2Z)id. 

Here, we take a simplified approach that shows how Dq of one-dimensional cellular fiame fronts can be 

related to e. By doing so, we neglect some subtle, but important, complications that arise in the statistical 
treatment of the Sivashinsky equation which are accounted for in the more sophisticated approach of 
Blinnikov and Sasorov (1996). In our simple model, we describe the cellular front as a hierarchy of cells on 
all length scales. If consecutive cell-generations are widely separated, we can assume that the thin flame 
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approximation is valid on each scale. The effective propagation speed Ui on scale li is then related to scale 
k^i < k by 

Ui = + eui_i . (4) 
Taking the continuum limit and integrating yields 

u, = uoe" . (5) 

We now assume that cell splitting occurs after a dilation interval S, so that k — SU^i. If the smallest 
unstable length scale is of the order of the Markstein length, /q = ^crit, we find U = S'^^crit and we can 
express the effective flame velocity in terms of the length scale I: 

^ ^ e/lnS 



Wccn(0 « Wlam ; • (6) 

V 'crit / 

If the flame speed is interpreted in a geometrical way, i.e., u(l) oc Ai, where Ai denotes the increased surface 
of a cellular front observed at the scale I, it follows that the surface area behaves like a fractal (Mandelbrot 
1983) 



Al = Alam [^] (7) 
\ 'crit / 

with the fractal dimension 

^ccu - 1 + 1^ • (8) 

Inserting (||) into (^) we find, in agreement with Blinnikov and Sasorov (1996), that the fractal excess 
of (^ is proportional to 7^ = (1 — /i)^. Specifically, we find Do = 7r^/241nS', and agreement with the 
authors' numerical results yields a dilation interval, S" ~ 4. 



2.4. Nuclear burning in the flamelet regime 

In the case of a cellular flame front driven purely by the LD instability, i.e. in the absence of gravity, 
there is no known upper limit for the largest scale of cell formation. However, if gravity is turned on 
the cell structure is no longer scale independent, which leads to the break-down of nonlinear stabilization 
(Khokhlov 1995). As soon as the process of cell disruption and bubble formation occurs on the largest 
scales, a cascade of turbulent velocity fluctuations is produced that continues on scales below the actual 
large scale flame instability (Niemeyer & Hillebrandt 1995b). After the establishment of the turbulent 
cascade, there exists a range of scales where burning is dominated by isotropic, fully developed turbulence. 
We will restrict our discussion to the conservative assumption that the production of turbulence is provided 
purely by large scale fluid instabilities, so that turbulent burning can be called "passive". Some thoughts on 
"active" turbulent combustion, where thermal expansion within the burning region is assumed to influence 
the properties of turbulence, are given in section ( [3.1.2 ). 

Cell formation ceases on the scale I where the turbulent velocity fluctuations v{l) become comparable 
with Uccii(0- Here, the time of front interaction with turbulent eddies becomes comparable to the eddy 
turnover time and, consequently, perturbations caused by turbulence grow to amplitudes comparable 
with the cell amplitudes. The cellular flame front is thus unstable on these scales (Zeldovich et al. 1980). 
Furthermore, numerical simulations of curved flames subject to shear show the breakdown of the nonlinear 
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stabilization mechanism if Wshoar ~ Ucdi (Niemeyer & Hillebrandt 1996). The transition between cellular 
and turbulent burning regimes is therefore marked by the Gibson scale, defined as (Peters 1988): 

W(^gibs) = "ccll(^gibs) , (9) 

Note that, due to the arguments above, we take eq. (|^) as an equality, which differs from earlier results 
(Khokhlov 1995) where the smallest turbulent scale was found to be « 500^gibs- Using the Kolmogorov 
scaling law v{l) oc l^^^ it follows that Zgibs scales with the third power of Uceii/v{L), where v{L) is, for 
instance, the magnitude of the turbulent velocity on the largest turbulent scale L. If we assume that near 
the beginning of the explosion this ratio is close to unity, implying that ^gibs ~ L, and later decreases owing 
to the decreasing flame speed and the build-up of shear in the RT mixing region, we find that the Gibson 
length decreases continually during the explosion. An upper bound for the intensity of turbulence on large 
scales is given by the "freezing out" of turbulent motions as a result of the overall expansion of the star 
(Khokhlov 1995). Using v{L) < 10^ cm s~^ at i « 10® cm and Uccii ~ uiam > 10^ cm s^^ (Timmes & 
Woosley 1992), we estimate that /gibs ^ 1 cm. Consequently, /gibs ^ ^th for all densities p ^ 3 x 10^ g 
cm~'^. For this reason, turbulent nuclear flame fronts in white dwarf matter at these densities burn in the 
"corrugated flamelet regime" (e.g., Clavin 1994). At lower densities, the flame front enters the so-called 



"distributed regime" that will be described in sections (2.6) and (3.2.5) 



The flamelet regime is characterized by laminar flame propagation on microscopic scales, while burning 
is determined purely by turbulence on large scales and therefore independent of microphysics. Fuel 
"digestion" occurs in an extended region behind the boundary of fuel and ashes, the so-called "turbulent 
flame brush" . The size of this region and the range of turbulent length scales adapt in a way that provides 
a flame brush propagation velocity that is decoupled from the laminar flame speed. On scales obeying 
/ ^ /gibs; one can express the effective flame speed wtur(0 and turbulent front width, dtnr{l) in terms of 
I and the turbulent velocity v{l) (Kerstein 1988, Clavin 1990), since both transport of burning fluid into 
the fresh material and fuel consumption inside the flame brush are limited by the eddy turnover time 
Tt{l) ~ OC l"^/^. Consequently, the relations 

Wtur(0««(0 , dtur(0«^ (10) 

are reasonable order of magnitude approximations that have been employed extensively in combustion 
research in the limit v{l) wiam- 



2.5. The Ray leigh- Taylor regime 

On the largest scales (/ > 10® cm) the flame dynamics is dominated by buoyancy of the hot, burned 
material surrounded by denser carbon and oxygen. The linearized problem is expressed by the RT instability 
(Chandrasekhar 1961). After a short period of exponential growth, the perturbation amplitudes become 
comparable to their wavelengths and the structure enters the nonlinear stage where interactions among the 
structures can occur, giving rise to merging and fragmentation of bubbles. Another important process is the 
appearance of KH unstable regions along the bubble surfaces. These produce turbulent eddies that spread 
out bubble tips and walls. Finally, the long time evolution shows continuing merging and fragmentation of 
rising bubbles, creating an increasingly turbulent mixing layer (Snider & Andrews 1994). 

Experiments of gas bubbles rising in vertical tubes flUed with fluid were performed to measure the 
asymptotic velocity of a single bubble (Davies & Taylor 1949). Best flts to the measurements were obtained 
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by the relation 

Vrt = B y/gcffl , (11) 
where B G [0.466,0.490] is a constant, ^off = Atg is the effective gravitational acceleration 
{At — (pu — Ph)/{Pu + Ph) ~ 0.5(1 — /i) is the Atwood number), and I, the radius of the tube. 
Layzer (1955) solved the problem analytically for a spherically symmetric tube and derived B — 0.511. 
Khokhlov (1995) derived a similar prescription for the propagation speed of RT unstable flames in open 
boxes from numerical simulations. While this result is true for a single length scale, we need to consider 
a range of RT-unstable scales, where the front creates a large number of bubbles with various radii. The 
so-called "Sharp- Wheeler Model" statistically describes the nonlinear stage of a multiscale RT-front (Sharp 
1984). It consists of a one-dimensional string of bubbles that are described by their radius, r^, and height, 
Zi. According to (|ll|), the heights grow as Zi — Wj-t. Neighboring bubbles merge if the difference between 
their heights exceeds the radius of the smaller bubble and thus the average bubble radius increases with 
time. Consequently, the average rise velocity from eq. (|l^) also grows. Numerical simulations of the 
Sharp- Wheeler-model (Glimm & Li 1988) show that the front asymptotically attains a constant acceleration 
that is proportional to g^s, so that the RT-mixing region grows at a rate of approximately 

few « 0.05geffi' . (12) 

This result is in agreement with full hydrodynamical simulations (Young 1984) and experiments (Read 
1984). Equivalently, we can say that the front advances into the cold material with a speed of 

Usw « O.l^offi ■ (13) 

Comparing (pT|), ( p^ ) and (^3|) shows that the maximum bubble radius evolves linearly with its displacement 
from the stellar center (i.e., the inner boundary of the RT mixing zone), /max oc r. 

In the context of supernova modeling, we are mainly interested in the burning velocity of the RT 



unstable flame brush. As in section (2.4), we can argue that the overall burning rate is limited by the 
fastest transport mechanism that mixes ashes and fuel. Invoking the arguments of the previous section, fuel 
consumption automatically adjusts to the speed of fuel contamination by burning blobs in order to provide 
a burning rate that is independent of microphysics. The highest speed for each single length scale I is now 
given by (pT|), yielding Uy-t oc Z^/^. A geometrical interpretation (eq. 0) of the flame speed in the nonlinear 
RT-regime therefore yields a fractal dimension of the flame surface, D^t — 2.5. In addition, we need the 
evolution the maximum bubble size, Zmax, as a function of time or radial displacement, like the one provided 
by the Sharp- Wheeler model. If we assume that bubble growth is purely governed by merging, and that 
burning is completed somewhere within the RT mixing region (this need not necessarily be the case), the 
effective burning velocity is given by the Sharp- Wheeler speed ([l3|). One dimensional supernova models 
using eq. ( p^ will be presented in section ( ^.l.l| ). 

Finally, we define the boundary between the turbulent burning regime dominated by Kolmogorov 
scaling and the buoyancy driven RT-regime by looking at the minimum time scale for self-interaction of 
flame structures with the size /, roughly given by Ts\{1) w V"rt(0 stated above, expansion of 

the star inhibits the growth of the largest structures into the fully nonlinear regime where they become 
isotropically turbulent. Expansion is characterized by the hydrodynamical time scale Tdyn ~ 0.1 s, so 
that self-interaction resulting in fully developed turbulence occurs on scales below Ztur/rt defined by 
Tsi(/tui/it) = Tdyn- Inserting equation (|l^) yields 

knr/rt = T^y„ ^off ~ 10^ Cm , (14) 

where we have used ^cff ~ 5 x 10^ cm s^^ as a typical value for the effective gravitational acceleration. 



- 8 - 



2.6. Summary of instabilities and their effects 

Beginning at the smallest dynamically relevant scale, the thermal flame thickness, Ith, the flame 
propagates with the laminar flame speed, uiami until the smallest cells appear at Zcrit ~ 100 Zth- Therefore, 
u{l) = const(/). In the cellular regime, u{l) oc /-^^oii-a ~ ^o.i^ Cellular stabilization fails when turbulent 
velocities become comparable with the effective cellular flame speed at Igihs- As a result of the Kolmogorov 
scaling law for turbulent velocity fluctuations and the assumption that the turbulent flame speed is 
determined by the eddy turnover time on every scale, we find that oc l^^^ in the fully turbulent 
regime. Above /tur/rt: the largest upward velocity on scale / is determined by buoyancy. Hence, u{l) oc l^^"^ 
according to (pj]). Fig. (|l]) is a summary of the scale dependence of the burning speed u{l). The piecewise 
scale- invariant burning regimes are represented by straight lines with different slopes in the log- log-plot. 
Three separate curves are shown, corresponding to the characteristic velocities and length scales at early, 
central, and late times of the deflagration phase (as represented by three decreasing density values). All 
graphs are based on the assumption that all the relevant instabilities and the transitions between them 
have reached a statistical equilibrium state. This is true after approximately one growth time of the largest 
considered scale (« 0.2 s). 

The curves in Fig. (0) are constructed as follows: First, the largest fully turbulent length scale km/i-t is 
computed from equation ( |l4|) with values for g taken from Khokhlov (1993a), while the burning speed at 



this scale is determined from equation (11). This point in the u-/-plane serves as the origin for the purely 
buoyant part with u{l) oc /^/^ extending to larger I, and for the Kolmogorov part reaching downward. The 
laminar flame speed is used as the second absolute point of each graph. Its value, as well as the thermal 
flame thickness, Ith, and the expansion factor, 7 = Ap/p, are taken from Timmes & Woosley (1992). At a 
length of /crit ~ 100 hh, the LD-instability and the subsequent formation of cells mark the transition to the 



scaling u{l) oc l^"-"-^ with L»coU ~ 2(1 + 0.87^) (Bhnnikov & Sasorov 1996; HJ). This fine is extended until 
it intersects the Kolmogorov line coming from above, which defines the Gibson scale Zgibs- Notice that, by 
virtue of this construction procedure refiecting our current understanding of flame dynamics in scale space, 
all burning properties above Zgibs are uniquely flxed by large scale phenomena and therefore independent of 
microphysics. 

The early phase of the explosion, where p « 10^ g cm~'^, is characterized by a high laminar flame 
speed and small thermal expansion, resulting in a very shallow slope of u{l) in the cellular regime. On 
scales of approximately 10 km, however, both effects are already dominated by the RT-instability and 
turbulent burning. The turbulent regime becomes more pronounced as wiam declines and the Gibson scale 
decreases {p ~ 10^ g cm~^, second graph). Meanwhile, the cellular part of u{l) becomes steeper owing 
to the increasing thermal expansion 7. In the simplified construction method outlined above, the highest 
turbulent velocities are assumed to be coupled directly to the RT- velocity at the transition length /tur/rt- 
Consequently, the model implied here "forgets" about all turbulence that has been built up by earlier 
RT- fluctuations, in contrast with the more realistic expectation that the level of turbulence decreases 
according to the almost negligible microscopic dissipation into heat. This assumption of instantaneous 
adjustment results in decreasing RT- and turbulent flame speeds as the density - and therefore g^e in 
equation (|lj) - falls. A phenomenological approach to account for the memory effect of accumulated 
turbulence in numerical simulations of SN la's has been proposed by Niemeyer & Hillebrandt (1995b). 

Fig. (^ also shows that our argumentation breaks down at the lowest density (p w 10^ g cm~'^), 
represented in the bottom graph. Here, the Gibson length, Zgibs is smaller than the thermal flame thickness, 
Ith, by almost four orders of magnitude. Therefore, the conditions for the flamelet regime (section 2A) are 
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Length Scale 1 [cm] 

Fig. 1. — Effective burning speed as a function of lengtfi scale for the early, main, and late phase of the 
deflagration. The shaded region represents the typical resolution of multidimensional simulations, A » 10 
km. 



- 10 - 




Fig. 2. — Regions of flamelet and distributed burning depending on fuel density and turbulent velocities at 

i pa 10 km. The lines depict Zth = ^gibs for different compositions ( Xc = 0.2, ■ • -Xq = 1, — Xq = 0.5). 

The shaded region shows the range of turbulent fluctuations induced by the RT-instability on large scales. 
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no longer satisfied, implying that laminar flame propagation ceases to occur on all scales (strictly speaking, 
this also implies that Zgibs is no longer a well defined quantity). Instead, turbulent eddies disrupt the 
burning region and dominate over conductive transport even on microscopic scales. This burning regime 
is sometimes termed "distributed" or "stirred combustion" (e.g., Pope 1987). We emphasize that all of 
the preceding discussion is based sensitively upon the assumption of microscopic flamelets and therefore 
becomes largely invalid after the transition to the distributed burning regime. Regions of distributed and 
flamelet regimes are displayed in Fig. (^) as a function of fuel density and macroscopic turbulent velocity 
i;(10km). The lines depict the transition defined by Ith — ^gibs for three different compositions. Decreasing 
carbon mass fraction is related to an earlier transition to distributed burning, since the flame thickness 
grows while the flame speed decreases. 



3. HOW THE WHITE DWARF EXPLODES 
3.1. Deflagration 

3.1.1. A simple deflagration? 

The simplest outcome to the exploding white dwarf problem would be a supernova in which the flame 
remained at all times and places subsonic and burned a sufficient fraction of the stellar mass to explode 
the star violently on the first attempt with no intervening pulsation. This is the basis for the common 
"deflagration" or "convective deflagration" model of which W7 (Nomoto, Thielemann, and Yokoi 1984) is 
a popular example. Numerous studies of this model (and others like it, e.g., Woosley, Axelrod, & Weaver 
1984; Woosley & Weaver 1994) give at least moderately good agreement with the observed light curve, 
spectrum, and nucleosynthetic requirements of common Type la events. The model does so by virtue of 
having approximately the correct proportions of intermediate mass elements (Si-Ca) to ^^Ni and the right 
final ^^Ni mass for a Chandrasekhar mass starting point. This result is, in turn, a consequence of a fiame 
speed that has been crafted to yield the desired result by a particular choice of transport algorithm (mixing 
length convection) and scale parameter (the mixing length as a function of time). Knowing approximately 
the desired result is useful, but can the flame speed in a real deflagration behave in the required way? 

Recent multi-dimensional calculations differ in their conclusions regarding this important issue. 
Khokhlov (1995) and Arnett & Livne (1994a) find that the flame moves too slowly to explode the star on 
the first try. Niemeyer & Hillebrandt (1995b) and Niemeyer et al. (1996) find that a prompt explosion is 
possible, albeit a weak one that makes little ^^Ni. A key difference in these calculations is the treatment 
of turbulence. All assume that the RT-instability is chiefly responsible for the production of turbulence, 
but while the former models assume that RT structures are isotropic and therefore implicitly employ an 
instantaneous adjustment of the turbulent cascade to the RT fluctuations, the latter simulations allow for 
a time delay between the production and dissipation of turbulent kinetic energy. They therefore include a 
"memory" for the history of the turbulent energy of each fluid element which, in some cases, gives rise to a 
build-up of turbulence above its local equilibrium. 

In the next sections, two physical ways of boosting the deflagration efficiency are discussed. A third 
possibility is simply that calculations having the necessary three dimensional resolution and low numerical 
viscosity have yet to be done. Were it not for drag (admittedly an unreasonable omission), a buoyant 
bubble floating in response to its density contrast (~40% typically when the flame is at 1000 km) and local 
gravity (~ 10^" cm s~^) would quickly achieve the sound speed. Assuming that the bubble dimension scaled 
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linearly with distance from the center of the star, a reasonable assumption in the bubble cascade model of 
Sharp and Wheeler (Sharp 1984; |2.5[), specifically taking the bubble radius, /,nax ~ 0.5r with r the distance 
to the stellar center and using eq. ([ll|), one obtains effective flame speeds ^1500 km s~^. This prescription 
implicitly includes drag and the result is consistent with what is seen by Garcia-Senz & Woosley (1995) in 
analytic calculations and by Niemeyer et al. (1996) in two-dimensional numerical models. A similar result 
can be obtained by directly using the Sharp- Wheeler burning velocity ( p^ ) which, for times 0.5 to 1 sec, 
is again of order 1000 to 2000 km s~^. Provided turbulence is capable of burning out or "digesting" all 
material internal to the leading bubbles, this speed would be enough to generate at least a mild explosion 
(consistent with Niemeyer et al. 1996). 

A simple one dimensional calculation illustrating this possibility is given in Fig. (^). A Chandrasekhar 
Mass white dwarf was ignited at its center and given an effective flame speed that, after some time (0.5 
s) during which the flame was assumed to move at a slow speed (approximately laminar), was assumed to 
be given by eq. (^). 5cff was calculated each time step at the flame boundary. Calculations were done 
using the one-dimensional implicit hydrodynamics code, Kepler (Weaver, Zimmerman, & Woosley 1978), 



also described in (3.2.1). This prescription assumes that the flame moves with a radial speed given for an 



ensemble of RT-unstable bubbles (Sharp 1984; 2.5) and that turbulent processes behind the leading edge 



can '"digest" all unburned fuel. The explosion is a healthy one, though short of intermediate mass elements. 

If for some reason the speed was greater, a stronger deflagration would result. Effects that have been 
left out of the simple Sharp- Wheeler speed (eq. ^ are bubble growth by burning, which causes even small 
bubbles to merge, the negative heat capacity of a mixture of ^^Ni and a-particles (Garcia-Senz & Woosley 
1995), and gradients in density and pressure. 



3.1.2. Active turbulent combustion? 



In section (2.4), we have neglected any feedback of burning itself on the intensity and spectrum of 
turbulence. In particular, we assume that the overall propagation speed of the turbulent flame brush is 
limited by the velocity of large scale RT-bubbles (equation |ll|). This need not be true. Owing to the 
limited range of observable scales, neither simulation nor experiment presently rule out the possibility that 
thermal expansion of burned material inside the flame brush significantly increases the strength of turbulent 
fluctuations on all turbulent scales. The energy available from nuclear burning is ~10^^ erg g~^, most of 
which goes into internal energy and bulk expansion of the star. Only a tiny fraction of this energy, ~10^^ 
- 10^*^ erg g^^, is injected into turbulence by the RT instability. Thus, the efficiency of any additional 
feedback mechanism that converts nuclear energy to turbulence need not be large. We therefore briefly 
discuss a possible mechanism for the feedback of burning on turbulence and its consequences for supernova 
explosions. 



Thermal expansion behind the laminar flame is the main reason for the LD-instability (section |2.2D . 
The growth of LD-perturbations in the absence of external fluctuations is stopped by cusp formation leading 
to a stable cellular flame structure. It is possible, however, that once the cellular structure is disrupted by 
the turbulent cascade from above, the production of specific volume contributes again to the growth of 
velocity fluctuations. Although local intersections of flame segments resembling cellular cusps still occur 
in the turbulent regime, their random orientation inhibits the formation of any self-stabilizing structure. 
Instead, it is likely that neighboring expanding regions accelerate the medium in some stochastic direction, 
thereby enhancing the intensity of velocity fluctuations. 
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Fig. 3. — Final velocity and composition in a simple one dimensional calculation of a deflagration in which 
the flame was constrained to move at a speed (in the comoving frame) of 0.1geg{t — to) with to = 0.5 s. The 
velocity is sampled at t = 10 s when it has taken on its coasting value. The final kinetic energy is 8.5 x 10^° 
erg. The mass of ^^Ni is 0.57 M© and there are 0.24 of other iron group isotopes. 
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In a first approximation, we assume that this effect is most efficient on the scale of the expanding 
regions themselves and therefore does not significantly couple turbulent velocities on different scales. As 
a consequence, only the fluctuation amplitudes would be affected, not their spectrum. The simplest case 
would then be a constant and scale independent growth rate, (j^tun of turbulent velocity fluctuations on 
each scale that depended only upon the expansion parameter 7, 

^-coMlXl), (15) 

giving rise to exponential growth: 

= v(Z,0)e'^'"('^)* . (16) 

Since the most plausible functional dependence of Wtur on 7 would be monotonously growing, this effect 
could readily account for the required delay of flame acceleration for a delayed detonation. Furthermore, 
the argumentation of section ( ^.4[ ) concerning the transport and burning time scales is still true, so that 
equation (^0|) still holds on every scale above the Gibson length (which now, of course, decreases according 
to the increasing strength of turbulence). More complicated, scale dependent feedback mechanisms would 
also alter the spectrum of turbulent velocities, giving rise to a fractal dimension of the flame surface different 
from D ^ ^3. 

A recent proposal by Kerstein (1996) gives strong support to the idea of active turbulent combustion. 
The author states that power-law growth of the turbulent burning velocity (as opposed to exponential 
growth expressed by eq. |l^) naturally occurs in the absence of a stabilizing mechanism, either by 
self-acceleration or through a statistical effect. However, the growth exponent is sensitive to the underlying 
model of flame dynamics and thus remains undetermined. 

As there is no natural limit for this growth in the incompressible regime, the velocity of turbulent 
eddies on large scales would eventually exceed that of rising RT-bubbles (equation |Tll). From there on, the 
properties of the flame brush would be independent of the RT-instability. As the turbulent velocities become 
close to sonic, compressibility effects grow more important and may give rise to a delayed detonation (see 



section 3.2.4). It must be noted that this self-enhancement effect can only occur after the cellular structure 
has been disrupted above the Gibson scale by the turbulent cascade from larger scales. This, together with 
the increasing gas expansion as the star expands, may explain the delay in the acceleration of the turbulent 
burning front necessary to account for SN la nucleosynthesis. We emphasize again that the arguments 
given here are speculation and need to be confirmed experimentally or numerically. Three dimensional 
simulations of burning in a pre-turbulizcd medium covering as many scales as reasonable are needed. 



3.1.3. Initial conditions - multi-point ignition 

As pointed out by Garcia-Senz and Woosley (1995) and demonstrated by Niemeyer et al. (1996), the 
efficiency of a white dwarf deflagration is sensitive to the manner in which the runaway is initiated. Multiple 
point ignition a few hundred km out from the center gives more powerful explosions than one ignited at the 
stellar center. This is because burning bubbles born off-center are more buoyant at birth due to the higher 
gravity and can travel farther before the star begins to disrupt. 

It is presently uncertain how far this trend can be continued. Obviously a white dwarf ignited 
simultaneously at at a large number of points scattered uniformly throughout its mass would have no 
difficulty exploding, but that would be artificial. The condition of simultaneity requires synchronization 
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of the burning in regions, that at the time they run away, are out of communication. The calculations of 
Garcia-Senz & Woosley indicate that it would take an extreme fine tuning of the temperature in order to get 
a convective bubble in the phase immediately preceding runaway to float more than ~200 km. Still a large 
number of 10 km bubbles, say, could be crowded on a sphere of this radius. If these floated a large distance 
after burning before acquiring sufficient velocity to drive lateral burning by the Kelvin Helmholtz instability 
and turbulence, the effective burning rate could become quite large (picture a dandelion gone to seed as the 
location of burned ash at this point). The time-dependent competition between lateral and radial burning 
is thus important. Three dimensional calculations to explore this possibility would be interesting. 

3.2. Detonation 

The simplest form of explosion, in so far as the hydrodynamics is concerned, would be detonation. 
This was the solution originally proposed by Arnett (1969), but numerous variations have been proposed 
since then that differ in where and when the transition from subsonic to supersonic burning occurs. In 
order to make this transition it is necessary that some volume sufficiently large to maintain a detonation 
burn in a time short compared to that required for sound crossing. Before discussing specific models, we 
begin by deriving a grid of critical masses capable of igniting and sustaining a detonation. We concentrate 
on densities between 10^ and 10* g cm~^ because this is the approximate range required for a transition to 
detonation that will provide ample intermediate mass elements. 

3.2.1. Critical masses for detonation 

The general idea of a critical mass for detonation follows Blinnikov & Khokhlov (1986, 1987), Khokhlov 
(1991a), and Woosley (1990); see also He & Clavin (1994) for a recent analysis of critical conditions for 
detonations. A propagating detonation wave must exert sufficient over-pressure to burn a region in a 
sound crossing time. If it does not, the shock will degenerate into a pressure wave and damp. The shock 
temperature for the self-consistent wave is given by the density and composition. This temperature gives a 
time scale for burning which when multiplied by the sound speed gives a rough estimate of the detonation 
wave thickness. Unless this distance is very small compared to the size of the region where the phase 
velocity of burning is supersonic, the detonation will damp. 

We have determined empirically a series of critical masses using the one-dimensional hydrodynamics 
code, Kepler (Weaver, Zimmerman, & Woosley 1978), which has nuclear physics and an equation of state 
appropriate for the problem. The nuclear network employed contained 19 isotopes from hydrogen to ^^Ni. 
Radiation transport was turned off and only the quadratic artificial viscosity term employed. Spheres were 
constructed of 100 zones with the mass of each zone smoothly increasing outwards (the first 30 zones 
had constant mass, the other 70 gradually increased logarithmically to a value 100 times larger than the 
central zone in zone 100). For the calculations at p < 10* g cm~'^, a linear (with respect to interior mass) 
temperature gradient was superimposed on the inner 28 zones with a central value of 3.2 xlO^ K falling 
to 4 X 10* K in zone 28. The rest of the zones also had a temperature of 4 x 10* K. This temperature 
gradient was such as to give a well resolved supersonic phase velocity for the burning in the inner zones for 
all spheres considered. The density was taken to be approximately constant in all zones and a boundary 
pressure applied to the sphere equal to that in the outer zones. For the p = 2 x 10^ g cm~^ case, the 
central temperature was assumed to be 2.8 x 10^ K declining smoothly over 55 zones to 7 x 10* K where it 
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remained constant for another 145 zones. A sample calculation is given in Fig. (Q). 

Burning was then turned on and the subsequent evolution followed for about 4000 times steps in each 
case (20,000 in the case of 2 x 10^ g cm~'^). Several dozen such spheres were modeled. A sampling of results 
is given in Table 1. In each case the critical mass is the lowest value for which detonation was achieved and 
maintained. The value given is the mass interior to zone 20, i.e., most of the material upon which the high 
burning temperature was imposed. In each case a calculation using 10 times less mass was also carried out 
and gave a failure. The length scale is thus resolved to about a factor of two. Because the results depend 
so sensitively upon the composition, it was not worthwhile to attempt greater resolution. 

For the case of 50% carbon and 50% oxygen our results agree well (within a factor of 10) with those 
of Arnett & Livne (1994b). However, we also explored the sensitivity to composition and found the mass 
fraction of carbon to be very important. At a density of 3 x 10^ g cm^"^, for example, the critical mass is 
5 orders of magnitude smaller for pure carbon than for a mixture of 50% C and O. The larger Q-value for 
the burning gives a higher burning temperature in the detonation and the reactions are very temperature 
sensitive. On the other hand, the mass is five orders of magnitude larger if the fuel is 40% carbon and 60% 
silicon, and for 35% carbon and 65% silicon the critical mass is larger than the star. We calculated mixtures 
of silicon and carbon because, later, we shall consider the possibility of igniting a detonation in a region 
where ash and fuel have mixed. At these densities, carbon and oxygen burn to intermediate mass elements 
of which silicon is more representative than nickel. Apparently, for each density, there is a critical mass of 
contamination by ash that makes detonation impossible. But if the ash remains a small component, the 
required critical masses are quite small. 

While the calculations were carried out for thoroughly mixed combinations, similar restrictions would 
apply to a region of thinly laminated, but unmixed fuel and ash. The critical quantity is the overpressure 
produced in a region large compared to the mixing scale. Since the ashes have long since yielded their 
overpressure to sustaining the subsonic expansion of the star, only the contribution of the fuel is important. 
If the mass fraction of fuel in the macroscopic region is low, the critical mass is large. 

These masses are for isolated hot spots surrounded by cold material. If the runaway occurs in an 
extended region that is itself close to burning on a sound crossing time, they could be smaller. 



We now discuss the possibilities for detonation in the approximate time order in which the detonation 
appears. 

First, there is the possibility of a prompt detonation (Arnett 1969). Critical masses for prompt 
detonation have been calculated by Blinnikov & Khokhlov (1986, 1987) and are very small (see also Table 
1 for p — 2 X 10^ g cm~'^), so there are an enormous number (^^10^^) of possibilities for detonation in the 
inner 100 km of the white dwarf that runs away. At issue, however, is the high degree of isothermality 
required to produce a supersonic phase velocity. 

Woosley(1990) derives the condition that the phase velocity of nuclear burning be sonic in the presence 
of a temperature gradient dT/dt: 



3.2.2. Prompt detonation 
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Fig. 4. — Microdetonation of a sphere of 50% carbon and 50% oxygen at 2 x 10^ g cm~'^. The detonation 
was intitated by a temperature gradient that gave a supersonic phase velocity in the inner 70 cm (10~^^ 
M0). The temperature fell smoothly in the initial model to a constant background of 7 x 10^ K at 87 cm 
(2.7 X 10^^^ M0). The density jump in the steady state shock was 1.5. The times sampled in the figure 
are 0.66, 1.9, 2.6, 3.2, 3.8, 4.3, 4.8, 5.4, and 6.5 xlO~® seconds. Velocity is the solid curve; pressure is the 
dashed curve. 
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Convection ceases to be efficient in transporting energy at about Tg = 0.7, when the convective cycle 
time equals the nuclear burning time, and this is a convenient point at which to evaluate the isothermal 
condition. The above equation gives a critical gradient, dT/dr ^ 0.001 K cm~^. Conduction will tend to 
smooth out temperature fluctuations on small scales. For p = 2 x 10^ g cm^'^ and Tg = 0.7, the conductive 
opacity is about 10~^ cm^ g^^ and the conductive length scale for 1 s, '-^lO m. The time scale to increase 
nuclear energy generation at Tg = 0.7 is ~10 s. Interestingly the critical mass for detonation (Table 1) is 
also a few meters in radius, so perhaps an approximately isothermal region of this size is not unreasonable. 

However, for detonation within a single critical mass of 10^^^ Mq (Table 1), the temperature could 
only vary from Tg = 0.7 by about one part in 10^" (if the conditions at Tg « 0.7 evolved without outside 
influence to runaway)! For larger length scales, up to 10*" cm, the restriction is less severe, but still very 
constraining, about one part in 10^. Further, this thermal condition must be maintained for an appreciable 
fraction of the nuclear time scale at Tg = 0.7, 

Tnuc = 0.15T9-20-2p-3.05 ^ 2O S (18) 

A prompt detonation becomes more likely if an only approximately isothermal region exists at higher 
temperature (eq. pT| ). For example, consider the conditions after the runaway has progressed and the 
hottest temperature anywhere in the core is Tg — 1.2. If this point is surrounded by say, 10 km of material 
in which the temperature is between Tg ~ 1.2 and 1.1, a detonation will occur. Similarly a point having Tg 
= 1.5 must be surrounded by 10"' cm of material between Tg = 1.4 and 1.5, etc. However, the nuclear time 
scale at these points are 5 x 10"'' s and 5 x 10"^ s respectively, and there is no efficient energy transport 
mechanism that operates on such short times. 

The actual thermal structure is not known with such precision. Small fluctuations will be amplified by 
the temperature sensitivity of the nuclear burning. Turbulence, energized by convection, has two effects. 
On the one hand, it acts to smooth out inhomogeneities. But, to the extent that the energy density in 
turbulence varies even slightly from place to place, different amounts of heat will be dissipated at the (very 
small) Kolmogorov scale. The convective velocity at Tg = 0.7 is approximately 30 km s~^. The energy 
density in this turbulence is about 10^^ - 10^^ erg g~^; the heat capacity of the gas, ~10^ erg g~^ K"^. 
The time scale for generating this energy is on the order of the convective turnover time 10 - 100 s. If 10% 
variations occurred in the dissipation rate in two regions, the temperature variation would be of order 10^ 
K. Though this does not prove that a prompt detonation is impossible, it suggests that it is much more 
difficult than the large number of critical masses implies. 

The chief evidence against prompt detonation is not so much the physics of igniting the detonation, 
as the fact that the models do not agree with observations. They produce too much ^"'Fe with respect to 
'^^Fe and no intermediate mass elements (e.g., Woosley 1990). The detonation wave may be unstable at 
high densities (p ^ 2 x 10'' g cm~'^; Khokhlov 1993b), but the instabilities predicted act on scales of order 
one cm and less. This would have no effect on the nucleosynthesis since even a laminar flame would burn 
through the residual fuel in less than an expansion time. If an instability were found that left behind regions 
of unburned material with the size of order 100 km, these regions would not be completely burned by 
conduction in the explosion. The model would resemble the deflagration ignited at many points throughout 
its volume discussed in the previous section. To our knowledge, no such instability has been proposed. 
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3.2.3. Atmospheric detonation 

The first models in which detonations were actuaUy calculated, rather than assumed (Nomoto, 
Thielemann, & Yokoi 1984; Woosley & Weaver 1986) occurred because a subsonic deflagration produced 
sufficiently strong pressure waves that their accumulation ahead of the flame led to compression adequate to 
ignite nuclear burning on a sonic time. This always occurred near the surface of the star in the steepening 
density gradient around 1.0 - 1.3 Mq. Such explosions made little intermediate mass elements, but did 
produce very high velocity iron synthesized as ^^Ni near the surface. The high velocities may be necessary 
to explain some supernova observations, especially for SN 1991T (Yamaoka et al. 1992), but the small 
synthesis of Si - Ca remains a problem for the typical SN la. The reason so little silicon is made in these 
models is that the flame must move an appreciable fraction of the sound speed in order to steepen into a 
shock in the outer layers. Such rapid burning consumes the white dwarf before it has had time to expand 
and is still at high density. 

To produce intermediate mass elements, it is helpful if the supernovae first burns relatively slowly and 
expands and then makes a transition to a detonation at a density ^10'' g cm~^. At such a low density 
detonation does not proceed to iron but to silicon. This defines the "delayed detonation" model of Khokhlov 
(1991abc) and Woosley & Weaver (1994). 



3.2.4- Delayed detonation 

First we consider the model of Khokhlov who attributes the transition to detonation to temperature 
fluctuations. These fluctuations arise from "non-uniform preheating of the gas" (or perhaps turbulence, 
though he emphasizes the former). Khokhlov points out that a spontaneous transition to detonation is 
frequently observed in terrestrial experiments (e.g., Lee & Moen 1980) involving turbulent flames. The 



underlying idea is similar to that discussed in (3.2.2). A single small region, but one larger than the critical 



mass, burns faster than sound because of an anomalously high temperature and shallow temperature 
gradient. Once a detonation is initiated, even a small one, it propagates through the rest of the star. 

However, such large localized temperature fluctuations are unlikely. So long as the flame retains a 



well defined surface (i.e., operates in the a laminar or "flamelet" regime; 2.4), high fuel temperatures exist 
only in the thin interface separating fuel and ash. Burning in this interface is very subsonic (defining, 
in fact, the laminar speed) and the mass in a flame thickness, much less than a critical mass. Previously 
existing, isolated temperature fluctuations in the carbon fuel will tend to be damped and quenched by the 
expansion. By very careful tuning, one might arrange to have a disconnected region run away at a late time 
(burning would occur at a rate that, for a long time, precisely balanced adiabatic losses as in Garcia-Senz & 
Woosley 1995). Additional tuning could give this region such a shallow temperature gradient that it burned 
supersonically. Not only does this seem unlikely, but the runaway would be much more likely to occur at 
an earlier time and there would be nothing special to delay the burning until the density declined below 
3 X 10^ g cm^'^, as is necessary if intermediate mass elements are to be produced. 

Temperature fluctuations might also be built up by the collision of sound waves or blobs of matter 
moving in random directions in the unburned carbon, but the energy density of these collisions is small. 
Typical velocity shears are 100 km s^^ implying an energy density of ~ lO^'' erg g~^. At a density of 3 x lO'' 
g cm~^, an internal energy addition of about 10^^ erg g~^ is required to raise the temperature above 2.0 
billion K where burning on a hydrodynamic time scale can occur. We conclude that a delayed detonation 
initiated by temperature fluctuations is unlikely. 
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A different sort of delayed detonation was discussed by Woosley & Weaver (1994; in lectures presented 
at Les Houches in 1990). This kind of model resembles more the atmospheric detonations of section (|3.2.3 ), 
but there is an attempt to "sculpt" the flame speed so as to provoke a detonation in a star with initially 
slow burning. In outcome the model is the same as that of Khokhlov, but the detonation occurs because of 
an accumulation of overpressure in a macroscopic region surrounding a topologically complex flame surface. 
Geometry plays the key role, not temperature fluctuations or gradients, and it makes sense to cast the 
model in terms of the fractal dimension of the flame. Some representative conditions illustrate the idea. At 
3 X 10^ g cm^'^, the critical mass for detonation is around 10^^"* Mq (50 meters for 50% carbon; Table 1). 
The turbulent flame brush (section ^.4| ) might be up to 200 km thick (i.e., it contains many critical masses; 
Niemeyer et al. 1996), the laminar speed of the flame, about 0.5 km s~^, the sound speed near 5000 km 
s~^, and the smallest unstable wavelength, about a centimeter. The 200 km thick region will then burn 
supersonically if the fractal dimension in that region exceeds 2.6. 

There are several difficulties here. First, the fractal dimension associated with a turbulent flame is 
well pegged to 2.3 - 2.36 based both upon simple scaling relations (Kerstein 1988, 1991) and experiment 
(e.g., Mantzaras et al. 1989, North & Santavicca 1990, Haslam & Ronney 1995). A value of 2.6 would 
require non-standard assumptions (though Mandelbrot, 1983, offers arguments that the fractal dimension 
of isoscalar surfaces in turbulence should be 2.5 to 2.66). It should be recognized however that D = ^3 is a 
statistical average. There may be regions having larger dimension and others with smaller ones. 

Why would this detonation occur preferentially at late times? At earlier times, the Gibson length is 
much larger and the structure of the flame brush correspondingly coarser. At 10^ g cm~^, for example, the 
laminar flame speed is about 50 km s^^. A flame brush of 200 km (though none would exist yet at this 
early time) would have a sound crossing time of about 0.02 s and the spacing between burning regions for 
supersonic effective burning would be 1 km. The Gibson length at these conditions is larger than 1 km; so 
too is the minimum unstable Rayleigh Taylor wavelength (Timmes & Woosley 1992). It would be hard to 
prepare such a layered structure without burning all the fuel in the process. Detonation would be difficult. 

A second condition required for a delayed detonation of this type is that the mass fraction of ash in the 
flame brush not be high. At 3 x 10^ g cm~'^. Table 1 suggests that detonation would not occur if the ash 
comprised more than 35% of the mass of the flame brush. 

We conclude that a delayed detonation of the kind proposed by Woosley & Weaver (1994) is possible, 
but its actual occurrence improbable. Whether the requisite fine structure can be set up is unknown. 

3.2.5. Detonation in the distributed regime 

Instead of mixing and wrinkling a thin flame sheet, turbulence, in the extreme, might mix heat from 
the ashes into cold unburned fuel. Rather than accumulate a critical surface area within a volume, the 
star might instead gradually build a critical temperature in a volume of fuel larger than in Table 1, as 
Khokhlov's model (1991abc) requires. This sounds straightforward, but can be difficult to arrange. Indeed, 
for densities above about 3 x 10^ g cm^'^, without pulsation, it is impossible. 

Above this density, the microphysics of the nuclear burning and electron conduction can be decoupled 
from all the instabilities we have discussed. Viewed on a sufficiently small scale (i.e., the conductive flame 
thickness), the flame is smooth and stable (Fig. |l|). Since the thickness of the flame is smaller than any 
critical mass in Table 1, no detonation can occur. 
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However, as the density drops below about 3 x 10^ g cm~'^ (Fig- H) nuclear burning proceeds at a 
slower rate and the flame becomes thicker. For a mixture of 50% carbon and 50% oxygen, Timmes & 
Woosley (1992) find a flame thickness of 0.5 cm at 5 x 10^ g cm~^ and 4 cm for 1 x 10^ g cm~^. At the same 
time, the scale of the smallest turbulent eddies that can turn over without burning (the Gibson scale) grows 
progressively smaller. For a turbulent energy density of 10^^ erg g^^ at 10® cm (Niemeyer & Hillebrandt 
1995b), the Gibson scale is 0.2 cm and 10"'* cm, respectively at 5 and 1 x 10*" g cm~^. Somewhere between 
these two densities a transition to a different kind of burning must occur. 

Under these conditions it no longer makes sense to speak of a flame propagated solely by electronic 
conduction. The burning is smeared out by turbulence; one enters the distributed regime of Fig. (|^). Heat 
can then be extracted from burning regions and transported to fuel. To make the transition to detonation 
a region of size larger than the critical mass (Table 1) must assume a temperature gradient shallower than 

with a peak temperature such that the nuclear burning time is much shorter than the stellar expansion 
time. We now illustrate, with specific conditions, how this might occur. 

First, a portion of the flame brush moves into the distributed region. We take a density 3 x 10^ g cm^"^ 
as representative. The turbulent flame brush is, at this point, a fine grained mixture of discrete phases - 
fuel and ash, fine grained, but still separated. Now, in places, the burning begins to go out as turbulence 
penetrates the flame sheet, homogenizes the composition, and reduces the temperature to some low mean 
value where burning is very slow. Recall the physics of the laminar flame (Timmes & Woosley 1992). There 
is a critical temperature in the flame where conduction balances nuclear energy generation. The nuclear 
time scale at this temperature times the laminar flame speed equals its thickness. As the flame enters the 
distributed region, turbulent eddies with a size comparable to the flame thickness (^ 1 cm) become fast 
enough to carry away heated fuel before it can burn. One might say that turbulence, for the first time, 
"enters" the reaction region and disperses it to larger length scales where, again, the burning time equals 
the (turbulent) transport time. For still lower densities and temperatures, and thus lower burning rates, the 
same level of turbulence is able to distribute the reaction zone on even larger scales. 

This leads to regions where the burning is temporarily quenched. As the overall density declines, these 
regions grow in size and, if the star becomes unbound, eventually encompass all the flame brush. These 
regions of suppressed burning are still coupled by turbulence, however. An important time scale to keep in 
mind is the characteristic turbulent time scale for a critical mass (Table 1). We shall assume, and justify 
later, that the fraction of ash in the mixture is small. Taking carbon equals oxygen equals 50% by mass, 
the critical mass is ^ 10^1" Mq, and its size 50 m. Its turbulent turnover time is then Tdit '^-^ 10 s (Fig. 
|l|). It is important that this is very much less than the hydrodynamical time, Tdyn = 446/pi/2 ~ 0.1 s, for 
the star to expand. 

If the mixture of fuel and ash burns faster than Tcrit , mixing will be incomplete and no detonation 
is possible. This requires that the temperature of the mixture be less than Tg — 2.2 (|l^), cooler if one 
considers larger regions containing many critical masses. On the other hand, the mixture must be able 
to resume rapid burning during the hydrodynamic time scale and that requires a temperature larger than 
Tg = 1.7. For intermediate temperatures of the mixture, a detonation is possible. The temperature in the 
ashes of a typical deflagration when the density is 3 x 10^ g cm^'^ is Tg = 4 to 5. The electronic heat 
capacity still dominates, so Cv oc T. Therefore mixing equal masses of fuel and ash gives a mixture with 
temperature Tash/2^^^, with little dependence on the temperature in the cold fuel. If we want the mixture 
to have a temperature of Tg ~ 2, for example, then each part of ash must be diluted with 5 parts of fuel if 
the unmixed ash has temperature Tg = 5 and 3 parts fuel if it has temperature Tg = 4. This justifies the 
use of the (ash-free) critical mass chosen above. The actual dilution factor must ultimately come from a 
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numerical study that we are unable to do right now. If it is too large or too small the model fails. 

A delay follows during which the temperature rises and the mixture runs away. For a detonation to 
occur, a region larger than the critical mass must have a nearly isothermal temperature distribution, (p^, 
about 20 K cm~^ if we continue to use our (very approximate) representative value of Tg = 2. This means 
is that once the mixture resumes burning, there can be no large region cooler than this peak value minus 
the temperature gradient times the size of the region (10^ K for 50 m). The region can, of course, be larger 
than the critical mass and the condition on the temperature fluctuations less stringent, but larger regions 
will have longer turbulent time scales and it will be harder to mix them without the fuel already burning. 



As in the case of prompt detonation ( 3.2.2| ), one has many opportunities for an improbable event, so 



the outcome is uncertain. However, because of the higher characteristic temperature (Tg « 2 rather that 
Tg « 0.7), the isothermal condition is not nearly so stringent, and with as many as ~ 10^^ critical masses in 
the flame brush, a detonation does not seem so unlikely. The congruence of the density where distributed 
burning can begin with that required for a detonation to make appreciable intermediate mass elements is 
also particularly encouraging for this model - another key difference with the prompt detonation model. 



3.2.6. Pulsational Detonation 

If inadequate fuel burns to disrupt the white dwarf on the first pulse, contraction will cause rekindled 
combustion and one or more pulsations will ensue. The burning conditions after each pulse will differ 
appreciably from the one before. If the pulse does not go below ^ 3 x 10^ g cm"'^, the boundary between 
ash and fuel remains intact, but the flame surface becomes dispersed throughout a greater fraction of the 
mass and more convoluted as well. If the pulse goes below 3 x 10 g cm~^, turbulence will mix both the 
composition and the heat across the interface. During the recompression, a large region of shared energy 
(or of dispersed flame) will reignite giving rise once more to the possibility of detonation, or at least very 
rapid combustion. 

Models of this sort have been studied by Khokhlov (1991b); Khokhlov, Miiller, and Hoflich (1993); 
Hoflich, Khokhlov, and Wheeler (1995) and, in two dimensions, by Arnett & Livne (1994ab) following early 
pioneering work on pulsational deflagration by Nomoto, Sugimoto, & Neo (1976) and Ivanova et al. (1974). 
More recently Woosley (1996) has demonstrated, using mixing length convection theory and a fractal flame 
in a one-dimensional model where pulsation was explicitly calculated, several possible outcomes including 



either detonation (the pulsational analog to 3.2.4 ) or an accelerated deflagration. However, no previous 
calculation has properly considered the critical role of turbulence. Without turbulence, the surface topology 
is not made much more complex during the pulse, nor is heat appreciably shared between fuel and ash. 
The heated region is instead confined to a narrow layer on the surface of a conductive flame. This is the 
definition of the flame thickness and, for a monotonically expanding supernova, this thickness is, at all 
times, thinner than the critical mass for detonation. It is possible, in principle, for a very large amplitude 
pulsation to lead to such a thick conductive flame that, during recontraction, the heated region encompasses 
a critical mass, an approach taken, for example, by Arnett & Livne (1994b). But for the same conditions, 
turbulence would dominate the heat transport. 

Observationally, the pulsational detonation model, or at least the single large pulse version, has a 
difficulty. If a large amplitude pulse is necessary for the initiation of the detonation, one would expect 
some white dwarfs to explode weakly without pulsing at all. Why should all initial pulses fall just short of 
producing an explosion? These supernovae with their faint broad light curves have, at the present time. 
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no observational counterpart. Perhaps they await discovery. Or perhaps, the explosion proceeds through a 
series of low energy pulses, only the last of which is always adequate to provoke detonation. In this latter 
case, the density 3 x 10^ g cni~^ is critical. Pulsations that do not go below this value would not lead to 
greatly accelerated burning because their flames would remain, on the small scale, laminar and dominated 
by conduction. The first pulsation to go below this density would experience appreciable mixing and heat 
sharing and greatly accelerated burning, perhaps detonation, on the next pulse. Burning on a hydrodynamic 
time scale would also be guaranteed since burning would be responsible for halting the recompression. 

What actually happens will not be known until realistic three dimensional models (with realistic 
ignition conditions; |3.1.3 ) have been calculated, both to show the failure of the first pulse and the degree of 



turbulent mixing during subsequent pulses. 



4. CONCLUSIONS 

Our first conclusion, perhaps not a very reassuring one to the observers, is that there are a lot of 
models for how a white dwarf explodes, none of which can be definitely excluded. However, our analysis 
suggests some models are more easily realized than others and calculations to clarify the situation. 

The simplest model, plain carbon defiagration, is an area where progress can and needs to be made. 
Currently one does not know, in the common case, whether the burning front, with all its instabilities, 
is able to unbind the star without an intervening pulse, which would exclude the pulsational detonation 
model, or not. If a radial front moves at the Sharp- Wheeler speed (|3|), the star becomes unbound, but 
this assumes efficient combustion by turbulence within the RT mixing layer The radial velocity 

of the bubble front is set by the non-linear Rayleigh- Taylor instability, which should advance at nearly 
the Sharp- Wheeler speed (0.1 gcst). Behind the leading edge of the bubble front, burning is enhanced by 
turbulence generated by the Kelvin-Helmholtz instability in shear fiows bounding the RT unstable blobs. 
It is not clear that fuel consumption proceeds at the same speed as the bubble front moves into the fuel. 
In other words, a stationary turbulent flame brush may never be established during the explosion, in which 
case the burning speed remains smaller than the Sharp- Wheeler speed. High resolution three dimensional 
calculations should ultimately clarify the issue. If the two speeds are comparable, an explosion seems likely, 
if only a mild one. The success of the first pulse is also very dependent upon how the runaway is initiated 
in the star ( p.l.l| , |3.1.3D . 



If the simple deflagration succeeds then, in order to agree with observations, it must have a higher 
speed than current multi-dimensional calculations suggest. We have discussed two possible ways of speeding 



it up - active turbulent combustion ( 3.1.2 ) that increases the turbulence intensity by thermal expansion 



within the flame brush, and extreme multi-point ignition (3.1.3). According to recent results for the scaling 
behavior of unconfined turbulent flames (Kerstein 1996), active turbulent combustion is a very promising 
field for future investigations, but its influence on the explosion mechanism remains speculative. 

Delayed detonations without an intervening pulsation are also not excluded, but those in the current 



literature have problems (3.2.4). The temperature fluctuations required to induce detonation by a 
"spontaneous burning" (Khokhlov 1991abc) are unlikely and the large fractal dimension and small minimum 
wavelength used by Woosley & Weaver (1994) are inconsistent with current views regarding turbulence. 
Indeed, in the current view, an effective burning speed, Utur, faster than the fastest turbulent motion, 
v{L), (which occurs on the largest scale L) is impossible (section ^^ ). This form of delayed detonation 
model survives for the time being, because the Kolmogorov mean field description of turbulence may 
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not fully describe our time dependent situation. The relation utm ~ '^{L) is an average. There may 
be appreciable fluctuations from this mean in isolated small regions larger than a critical mass (we are 
currently investigating this). Active turbulent combustion ( [3.1.2| ) might enhance the fuel consumption rate 
sufSciently. Another possibility is a burning geometry that favored only radial growth early on, but rapid 
non-radial combustion at late times. 

A more likely possibility, which we are suggesting here in the astrophysical context for the first time, is 
a transition to detonation as the star expands, its density declines, and it enters the regime of distributed 
burning (Fig. ||; |3.2.5 ). This kind of explosion has several appealing characteristics. First, the transition to 



detonation occurs as a deflagration is dying. This naturally leads to the desired pre-expansion of the star 
and rapid burning at late times. In particular, abundant intermediate mass elements can be synthesized at 
densities below the transition density, ~ 10^ g cm~^. If the detonation ignites, it will naturally give a very 
energetic explosion, currently a problem for the deflagrations. Most importantly it does these things using 
credible, definite physics, which, though uncertain, can be tested by numerical modelling. At this point 
its greatest uncertainties are the actual dilution factor for fuel and ash (a value of one part ash to several 
parts fuel is optimal) and whether the necessary isothermal conditions can be set up in the mixed medium 
at Tg « 2. For these reasons, this is the model that that we favor at the present time. The observational 
consequences of this form of detonation should be the same as for other forms of (parameterized) delayed 
detonation models. 

If the detonation fails to catch on the first try, and if simple deflagration fails, then pulsational 



detonation (3.2.6) offers an appealing alternative, but not one without its own puzzles. The deflagration 
must fail after having burned an appreciable mass so that the dwarf experiences a large amplitude 
oscillation. This requires some tuning so that the star does not become unbound and yet still expands 
enough to mix. Why do we see no supernovae that just barely exploded? Here again distributed burning 
may help. The detonation only lights once the pulsation has sufficient amplitude to go to 10^ g cm~'^ at 
which point efficient mixing occurs. The star could approch this limiting density by one or more pulsations. 

Because fundamental physics does not yet preclude several qualitatively different outcomes from very 
similar starting points, it is possible that they all happen to some extent. Even within a single class 
of model there is room for considerable variation owing to different ignition conditions, variable ratios 
of carbon to oxygen (which set the underlying conductive speed as well as affect the critical detonation 
mass), and the uncertain transition to detonation. Such diversity may be necessary to understand such 
distinctively different Type la supernovae as SN 1991bg and SN 1991T. 
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Table 1. Critical Masses 



p 


12c 


160 


28Si 


56Ni 


Mass 
[M©] 


Radius 


[109 K] 






0.5 


0.5 








10-10 


2 km 


3.5 


4 X 10-4 




0.5 





0.5 





10-8 


10 km 


3.3 


2 X 10-3 


3 X 10^ 


1.0 











10-19 


1 m 


5.5 


2 X 10-^ 




0.5 


0.5 








10-14 


50 m 


5.2 


1 X 10-5 




0.5 





0.5 





10-10 


1 km 


4.3 


2 X 10-4 




0.5 








0.5 


10-10 


1 km 


4.1 


2 X 10-4 




0.4 





0.6 





10-^ 


50 km 


3.9 


1 X 10-2 




0.25 


0.25 


0.5 





10-4 


100 km 


4.0 


2 X 10-1 




0.35 





0.65 





> 10-1 










1 











10-20 


40 cm 


6.5 


1 X 10-^ 




0.5 


0.5 








10-18 


2 m 


6.2 


4 X 10-^ 




0.5 





0.5 





10-12 


150 m 


5.5 


3 X 10-5 


2 X 10^ 


0.5 


0.5 








10-18 


70 cm 


9.7 


7 X 10-8 



^Sonic crossing time assuming a sound speed of 5000 km s i; (10,000 at 2 x 109) 
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